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ABSTRACT 


A  stochastic  stellar  atmosphere  model  which  simulates  the  behavior  of  pho- 
Ums  in  the  atmospheres  of  stars  by  following  individual  photons  is  described  and 
tested.  The  model  aocoimts  for  stellar  turbulence,  atmospheric  curvature,  stellar 
rotatkm,  oolHsional  excitation,  coUisional  de-excitation,  and  stellar  winds.  Testing 
shows  that  this  new  method  of  stellar  radiative  transfer  is  successful  in  simulating 
all  of  these  processes.  Particular  results  of  interest  include  simulation  of  solar  limb 
darkening  and  spectral  synthesis  of  three  different  solar  absorption  lines;  Ti 
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I.  INTRODUCTION 


Rec«nt  advADCM  in  •peclroacopy  hav«  l«<i  to  a  heii^tened  interest  in  radiative 
tranaier  the(»y.  As  R  P.  Kudhtski  stated  (Ref  1] 


If  o«ie  looJu  at  all  the  published  papers  resulting  from  extremely  expra- 
sive  space  expeiiments  or  extensive  ground  based  observations  us^ 
large  telescopes  then  a  <hsturfain|dy  large  fraction  contain  simply  line 
identifications,  radial  velocities,  and  only  qualitative  arguments.  In  view 
of  all  the  financial  and  engineering  eSori  that  goes  into  the  development 
of  our  modem  instruments  this  is.  of  course,  not  enough. 


What  is  tacking  is  a  comprehensive,  workable,  quantitative,  radiative  transfer 
theory  with  which  to  extract  as  much  information  as  possible  from  available  stellar 
spectra.  The  principle  problem  Csced  by  modem  radiative  transfer  theory  is  the 
treatment  (rf’  ooo*tocaJ  ibennodynamic  equilibrium  (NLTE)  portions  of  stellar  atmo¬ 
spheres.  A  few  examples  are:  B*  and  A-type  supergiant  stars  where  NLTE  theory 
is  needed  in  the  photosphere  and  the  srind  region,  solar-type  stars  where  strong 
absorption  lines  occmt  in  the  upper  jdiotospbere  and  the  chromosphere;  pre-main 
sequence  stars  with  hi|d>  mass  loss  and  disks  that  protrude  into  the  photosphere; 
and.  probably  the  moot  complex  case,  0-type  stars,  where  all  LTE  assiunptions 
breakdown  completely  even  deep  in  the  photosphere. 

Many  diflerent  themetical  techniques  have  been  applied  to  the  problem  of 
radiative  tnmsfer  with  varying  degrees  of  stMoess.  In  this  paper,  a  new  approach 
which  fidlows  individual  photons  throti^  the  stellar  atmosphere  of  the  star  is 
examined. 
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A.  BACKGROUND  ON  RADIATIVE  TRANSFER  THEORY 


Radiative  transfer  theory  is  the  study  of  how  photons  propagate  through  a 
medium.  The  theory  has  its  origins  in  astronomy  and  the  study  of  star  atmo¬ 
spheres  and  spectra.  But  there  are  many  additional  military  and  civilian  applica¬ 
tions  of  radiative  transfer  theory  (i.e.,  fusion  bomb  triggering,  nuclear  radiation 
transfer,  atmospheric  sciences). 

To  describe  the  interaction  of  matter  and  lig^t  it  is  necessary  to  characterize 
the  state  of  the  matter  as  well  as  that  of  the  radiation  field.  T^ically,  the 
assiimption  of  thermodjmamic  equilibrium  (TE)  or  local  thermodynamic  equilibrium 
(LTE)  is  made  to  simplify  this  characterization.  These  assumptions  significantly 
reduce  the  complexity  of  many  of  the  radiative  transfer  computations.  A  quick 
review  of  these  assiunptions  is  in  order.  TE  is  the  most  restrictive  case.  A  gas  in 
TE  will  satisfy  the  following  conditions: 

1.  Each  species  of  gas  particle  will  have  a  Maxwellian  velocity 
distribution  specified  only  by  its  mass  and  a  temperature; 

2.  The  distribution  of  the  excited  states  of  each  ion  will  be 
acciuutely  described  by  Boltzman’s  law  specified  only  by  the 
energies  of  excitation  and  a  temperatiire; 

3.  The  distribution  of  the  ionized  states  of  the  atoms  is  accur¬ 
ately  described  by  Saha’s  equation,  again  specified  only  by  the 
energies  of  ionization  and  a  temperature;  and 

4.  Radiation  intensity  is  accurately  described  by  Planck’s  function 
specified  only  by  a  temperature. 

It  is  important  to  note  that  in  TE,  the  same  temperature  variable  is  utilized 
in  all  four  relationships  given  above. 
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TE  can  only  be  achieved  in  an  iadated  eystem.  to  that  the  gae  can  not  inter¬ 
change  any  energy  or  mass  with  ite  eumnindinga.  It  ie  obvious  that  this  cooditioD 
is  rather  stringent  and  difficult  to  meet  in  a  stellar  atmosphere. 

Since  true  TE  cannot  be  fully  justified,  the  lees  strinfent  LTE  is  usually 
invoked.  In  this  second  case,  the  gas  wtU  have  TB  characteristios  over  a  local  tem¬ 
perature  zone.  TTiis  assximption  is  justified  when  collisions  among  atoms  of  the  gas 
occur  rapidly  enough  to  be  the  dominate  cause  of  esdtation.  When  this  occurs,  the 
kinetic  temperature  of  the  atoms  in  the  local  region  can  be  used  to  accurately 
describe  the  gases'  atomic  {Hoperties.  The  fourth  assumptioo  will  not  hold  for  the 
LTE  case  because  radiation  escaping  from  deeper  regions  of  the  star  cannot  be 
described  by  the  kinetic  temperature  of  the  local  region;  therefore,  the  radiatam 
field  can  deviate  arbitrarily  firom  the  Planck  distribution. 

The  LTE  assumption  is  a  good  deecriptioo  of  the  state  of  the  matter  deep  in 
the  photosphere  of  a  star,  where  the  migority  of  the  energy  is  transfom-d  via  con¬ 
duction  or  convection.  This  is  for  two  reasons:  1)  due  to  the  high  opacity,  little 
radiation  (energy)  can  escape  directly,  thus  preserving  energy  conservation;  and.  2) 
use  of  Boltzman's  etpiation  is  justified  because  collisions  are  the  main  mechanism 
for  atomic  ezdtaticm  thus  allowing  the  esdtation  to  be  described  by  a  single  tem¬ 
perature.  The  LTE  assumption  begins  to  break  down  as  the  oolhsiooal  transitioD 
rates  of  the  ions  approach  the  radiative  transitioo  .ates.  This  decrease  in  oolb- 
sional  rates  relative  to  the  radiative  rates  occurs  as  the  photons  work  their  way  up 
through  the  atmosphere  during  which  the  densi^  of  the  gas  and  the  kinetic  tem¬ 
perature  ai  tbe  atoms  decreases  as  energy  is  radiated  away. 

When  the  LTE  assumption  is  no  longer  valid,  the  state  of  the  gas  can  be  cal- 
nilated  by  statistical  equibbrium  (SE)  methods.  This  approach  balances  tbe 
nximber  of  transitions  dtu  to  all  proettmt  into  and  out  of  a  given  atomic  level  This 
method  is  very  difficult  given  the  large  ntimber  of  energy  levels  and  ionization 
states  available  to  an  atom.  Tbe  method  of  statistical  equilibrium  is  always  valid 
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•ad  Um  ••Aumptkxu  of  LTE  and  TE  ftr»  oaljr  tpedal  case*  (d‘S£. 

With  aa  uadarataadiof  ci  tha  charactahsation  of  tha  aiedium,  the  radiative 
traaafor  •<{uatioo  ia 


dtv 


■  f  *  "■ 


(1) 


whara;  /  «  >  iaitial  iateoaity  aa  a  fuactkm  of  firaqiiaacy, 

■  aourca  (uactioa  aa  a  fuaction  of  firaquaacy,  aad 
%,  m  optical  depth  aa  a  fuactum  of  firaquaacy. 

Xy,  aad  are  given  by 


L 


and 


(2) 


Sy  S 


h. 


(3) 


where:  dx  ■  geometrical  depth. 

L  a  path  length, 

p  a  deaaity, 

cm^ 

Ky  a  extinctkm  coefficient  ( - ),  and 

g 

jy  a  emission  coefficient  (  -  -  ). 

g  rad* 

Hie  simplicity  of  equation  (1)  is  deceiving.  The  complexity  of  the  radiative 
transfer  problem  lies  in  determining  the  properties  of  Ky  and  y'y  .  Ky  and  jy  are 
dependent  on  the  intensity  which  renders  equation  (1)  a  first  order  non-linear 
differential  equation  vastly  complicating  the  problem. 


Radiative  transfer  theory  has  been  successful  in  explaining  many  features  of 
stellar  spectra.  There  are  some  important  problems  in  which  the  restrictive 
assumptions  of  this  theory  fail  to  adequately  describe  the  stellar  atmosphere. 
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1.  NLTE 


Some  of  the  more  important  spectral  lines  (i.e.,  strong  absorption  lines 
due  to  H,  He,  and  metals)  have  non-LTE  (NLTE)  characteristics.  This  is  because 
the  line  centers  are  formed  high  in  the  atmosphere  where  the  density  and  kinetic 
temp>eratiu:e  are  low  and  radiative  processes  dominate  over  coUisional  processes. 
This  NLTE  behavior  is  not  restricted  to  exotic  high-temperatruie,  high-luminosity 
stars;  the  sun  contains  several  NLTE  lines  of  interest.  The  importance  of  under¬ 
standing  NLTE  line  formation  can  not  be  over  emphasized.  Spectral  measure¬ 
ments  of  faint  stars  cally  can  onlv  detect  strong  absorption  lines  over  back¬ 
ground  noise.  Many  of  these  absorption  lines  are  NLTE  in  natiire.  Therefore  an 
imderstanding  of  this  NLTE  line  formation  region  is  of  paramount  concern. 

2.  Winds 

Radiatively  driven  winds  in  stellar  atmospheres  are  an  additional  prob¬ 
lem  for  the  LTE  assumptions.  Many  stars,  particularly  high-temperature,  high- 
Ixxminosity,  very  young,  and  very  old  stars,  have  rapid  mass  loss  (typically  1  earth 
mass/yr)  via  radiatively-driven  winds.  These  winds  produce  expanding  stellar 
atmospheres  where  the  LTE  assumptions  quickly  break  down.  These  expanding 
atmospheres  also  preferentially  Doppler  shift  line  profiles,  which  complicates  the 
radiative  transfer  problem.  Loss  of  matter  must  also  be  considered  when  analyzing 
the  validity  of  the  LTE  assumption.  Boltzman’s  equation  requires  a  closed  system 
in  which  the  total  mass  of  the  system  is  conserved.  Rapid  mass  loss  is  in  violation 
of  this  requirement. 

B.  CURRENT  STELLAR  ATMOSPHERE  MODELING 

Stellar  atmosphere  modeling  is  an  attempt  to  reconstruct  the  structure  of  the 
stellar  atmosphere  based  on  spectral  observations  and  physical  laws.  The  complex¬ 
ities  of  the  dynamics  of  the  atmosphere  as  well  as  the  complexities  of  radiative 
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transfer  theory  prohibit  direct  measurement  of  many  of  the  physical  parameters. 
Therefore,  a  model  is  hypothesized  to  generate  the  observed  spectral  information. 
Once  the  model  has  been  constructed  which  reproduces  the  observed  information, 
additional  observations  can  be  made  to  further  test  the  model.  While  this  may  not 
produce  unique  atmospheric  parameters,  stellar  atmosphere  modeling  has  become 
an  indispensable  tool  in  the  analysis  of  stellar  spectra.  Most  models  contain  the 
following  simplifications: 

1.  LTE  is  assumed  throughout  the  atmosphere; 

2.  Plane-parallel  geometry  is  assumed  ignoring  atmospheric  ciirvature; 

3.  Hydrostatic  equilibrium  is  assumed  neglecting  any  mass  losses; 

4.  Fine  structure  is  negligible  (i.e.,  granulations,  flares,  siinspots);  and 

5.  Magnetic  fields  are  excluded. 

Stellar  modeling  can  be  divided  into  three  distinct  types:  continuum,  line, 
and  spectral  synthesis.  Each  method  will  be  examined  further  in  the  next  three 
sections. 


1.  Continuum  modeling 

Continuum  modeling  seeks  to  reproduce  the  characteristics  of  a  stars 
continuum  radiation.  The  source  fimction  is  calculated  firom  the  Planck  distribu¬ 
tion  at  a  given  depth  of  formation.  This  quantity  is  referred  to  as  the  surface  flux. 
The  extinction  term  from  the  radiative  transfer  equation  (eq  1)  is  based  on  contin¬ 
uum  absorption  coefficient  calculations  and  laboratory  measurements  of  atomic 
properties  (i.e.,  H  ,  H~  ,  electron  scattering).  The  standard  radiative  transfer 
integral  is  then  solved,  producing  the  calculated  spectrum.  The  success  of  these 
types  of  models  is  remarkable  considering  their  simplistic  assumptions.  Figure  1 
shows  a  fit  according  to  Melboure  [Ref.  2]. 
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Continuiiin  model  comparison  to  Observations 


Figure  1:  The  observed  energy  distribution  in  51  Peg(G5V)  corrected  for  line 
absorption  compared  with  model  calculations  [Ref.  2], 

Note  that  the  observed  data  has  been  adjusted  to  account  for  line  absorption  at  the 
shorter  wavelengths. 

The  continuum  modeling  technique  contains  three  mqjor  problem  areas: 
1)  absolute  continu\im  flux  calibration  (i.e.,  where  is  the  continuum  located),  2) 
NLTE  effects  on  the  temperature  gradient,  and  3)  interstellar  reddening.  The  first 
two  of  these  problems  is  intrinsic  to  the  continuum  modeling  technique  and  not 
easily  solved;  interstellar  reddening  complicates  comparisons  to  observed  contin¬ 
uum. 

2.  Line  modeling 

Line  modeling  efforts  have  been  directed  at  determining  the  shape  of  line 
absorption  profiles.  With  accurate  profiles  much  can  be  learned  about  stellar  atmo¬ 
spheres.  The  line  absorption  profile  is  dependent  on  many  important  atmospheric 
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parameter!  (i.e..  preeaure,  temperature,  radiatiao  6eld.  rotation,  winda,  turbulence, 
depth  of  line  formation,  etc.).  It  ia  the  interaetioo  of  theae  parameter!  that  {ho- 
duce!  the  abaorption  profile  and  thua  reveala  inlbnnatioo  about  the  condition!  at 
formation  of  the  line.  In  addition,  line  modeling  ivovidea  infurmation  about  the 
proceaaea  that  produce  variouB  line  profilea. 

The  radiative  tranafer  equation  (eq  1)  ia  the  baaia  of  Umae 
But  the  central  problem  remain!  how  to  <4Aain  an  accturate  aouroe  function.  Typi¬ 
cally,  Sy  ia  approximated  by  a  oombinatioo  of  a  line-aource  function  and  a 
continuum-aouroe  function.  The  Iine-!0\iroe  function  ia  calculated  from  LTE  am- 
siderationa  and  the  continuum-aource  function  ia  calculated  firom  aurface  flux  calcu- 
lationa.  Ibeae  calculation!  have  produced  good  reaulta  when  compared  to  actual 
absorption-line  profiles.  Figure  2  and  3  show  the  reaulta  of  line  mn^^ling  done  by 
Carbon  and  Gingerich  [Ref.  3],  and  Curtia  and  Jefleriea  (Ref  4)  respectively. 

Line  models  compared  to  observations 


Figure  2:  The  observed  Hy  profile  in  6  Leo  compared  with  LTE  line  profile  cal¬ 
culations.  The  model  calculations  fail  to  correctly  model  the  line  cores  [Ref. 
3]. 
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Figure  3:  The  obeerved  Na  D  line  profile  compared  with  LTE  line  profile  cal- 
culationa.  Again  the  model  calculations  fail  to  correctly  model  the  linft  cores 
[Ref.  41. 

The  depths  of  the  observed  profiles  are  deeper  than  those  calculated  by 
line-modeling  theory.  This  is  one  of  the  inadequacies  of  the  current  line-modeling 
theory.  Gray  dies  the  efTects  of  turbulence  and  NLTE  population  levels  as  typi¬ 
cally  deepening  line  cores  [Ref.  5).  Indusion  of  these  phenomena  improves  line  fits 
markedly  but  some  differences  still  remain.  The  line-modeling-technique  limita¬ 
tions  are  not  only  confined  to  NLTE  efTects  at  line  cores.  Weak-line  calmlations 
also  fail  to  match  theoretical  curves  of  growth.  Ciuwes  of  growth  describe  the 
theoretical  growth  of  an  absorption  line  as  a  function  of  abundance,  temperature, 
and  pressure.  Weak-line  calculations  fail  to  match  these  cxirves  although  indusion 
of  turbulence  improves  the  fits. 

3.  Spectra]  sjmtheals 

Spectral  synthesis  is  a  method  of  spectral  analysis  developed  to  calcu¬ 
late  the  chemical  abundance  in  stellar  atmospheres.  High-resolution  observations 
of  spectral  lines  are  compared  to  the  synthesized  lines.  Hie  individual  abundances 
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of  the  atoms  are  adjiisted  until  a  good  fit  in  the  synthesized  line  is  achieved,  thus 
determining  the  chemical  composition  of  the  star.  The  intensities  of  the  spectral 
lines  are  calculated  by  evaluating  the  standard  radiative  equation  (eq  1).  Tlie 
results  of  some  spectral  synthesis  are  shown  in  Figure  4  and  5. 

Spectral  synthesis  compared  with  observations 


Figure  4:  An  observed  portion  of  the  solar  spectrum  compared  with  Ross  and 
Aller  spectral  synthesis  [Ref  6]. 

Figure  4  is  a  comparison  of  computer  synthesized  lines  to  observed  solar  lines  [Ref 

6] .  The  fit  is  fairly  good,  although  turbulence  effects  have  been  added.  Figure  5  is 
the  synthesized  spectrum  of  a  carbon  star  (TX  Psc)  compared  to  observations  [Ref 

7] .  Although  the  fit  shows  significant  deviation  from  observations,  Figure  5  shows 
how  even  very  crowded  regions  of  the  spectrum  can  be  synthesized. 

This  third  approach  contains  many  of  the  same  problems  as  the  line¬ 
modeling  method.  The  emphasis  in  this  method  has  been  on  the  problem  of  line 
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Figure  5:  The  observed  spectrum  of  TX  Psc  compared  with  Synthesized  spec¬ 
trum  [Ref.  7]. 

blending  in  crowded  regions  and  not  on  the  line  profile  itself. 

4.  Limitatioiu 

Although  these  modeling  methods  have  been  very  successful,  each  does 
possess  some  limitations  as  noted.  In  addition,  all  of  these  methods  contain  four 
generic  problems:  1)  the  assumption  of  a  plane-parallel  geometry,  inappropriate  for 
large  extended  stellar  atmospheres;  2)  they  lack  an  easy  method  by  which  to  incor¬ 
porate  stellar  winds  and  rotation;  3)  Line  profiles  and  continuum  models  fail  to 
reproduce  this  observed  limb  darkening  behavior*  across  the  solar  disk;  and  4) 
LTE  atmospheres  are  assiuned. 


*  ObMTvationa  of  the  lun’e  disk  reveal  that  intensity  decreases  as  observations  near  the 
solar  limb. 


11 


5.  Other  approaches 

In  an  effort  to  solve  the  NLTE  problem,  Kudritzki  [Ref.  8]  used  the 
method  of  linearization  proposed  by  Auer  and  Mihalas  [Ref  9]  to  develop  a  com¬ 
puter  code  designed  to  model  the  structure  of  the  NLTE  atmosphere.  The  method 
accoimts  for  all  possible  coupling  between  variables  at  all  different  depths  within 
the  constraints  of  statistical,  hydrostatic,  and  radiative  equilibrium.  This  is  essen¬ 
tially  the  statistical  equilibrium  method.  To  show  the  complexities  of  these  calcula¬ 
tions,  the  O  II  Grotrian  diagram  used  by  Becker  and  Butler  [Ref  10, 11]  is  shown 
in  Figure  6. 

Grotrian  diagram 


Figure  6:  0  II  Grotrian  diagram  showing  the  200  different  transition  which 
must  be  calctilated  [Ref  10,  11]. 

The  computer  code  must  be  given,  or  calculate,  all  atomic  parameters  and  calculate 
all  transition  rates  into  and  out  of  these  200  different  atomic  energy  levels. 
Kudritzki  [Ref  12]  has  demonstrated  that  a  large  number  of  transitions  are  neces¬ 
sary  for  acceptable  accuracy.  With  accurate  atomic  parameters,  the  res\ilts  of  this 
method  can  by  highly  accurate  (Figure  7).  Discrepancies  in  the  line  core  were 
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corrected  for  by  Herrero  [Ref.  13]  by  including  new  coUisional  broadening  terms 
and  coupling  between  H  and  He  11. 


NLTE  line  profile  compared  to  observations 


•]6  0  S  16  -16  ~6  0  e  IS 


AX 

Figure  7:  Observed  Hy  profile  of  the  O  9.5  V  star  t  Sco  compared  to  NLTE 
calculations  [Ref.  9]  (left)  and  recent  improvements  by  Herrero  [Ref.  13) 
(right). 

In  an  effort  to  explain  multi-line  wind  transfer  problems,  Abbott  and 
Lucy  [Ref.  14]  developed  a  Monte  Carlo  technique  for  the  exploration  of  stellar 
wind  effects  on  the  spectra  firom  stars  with  strong  stellar  winds.  'Hie  atmospheric 
parameters  used  by  Abbott  and  Lucy  were  firom  Kudritzki,  Simon,  and  Hamann 
[Ref.  15].  Their  technique  follows  energy  packets,  which  may  be  as  small  as  /tv, 
firom  the  top  of  the  photosphere  to  their  emergence  through  the  stellar  wind  por¬ 
tion  of  the  upper  atmosphere.  In  this  way  wind  effects  can  be  simulated  and 
wind-driving  forces  explored.  Results  of  this  technique  are  shown  in  Figure  8  and 
reveals  general  agreement  with  observations  of  C  Pup  in  the  ultraviolet. 
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Wind  simiilatioD  compared  to  observations 


Figure  8:  Abbott  and  Lucy  fit  to  observed  ^  Pup  spectrum.  Prominent  features 
are  identified  on  the  top  axis  [Ref.  141. 


C.  STOCHASTIC  PHOTON-FOLLOWING  SD^IULATION 

1.  Motivation  for  SAM 

In  view  of  the  limitations  with  current  atmospheric  modeling  and  radia¬ 
tive  transfer  methods  a  new  approach  has  been  explored  by  W.  B.  Weaver  at  the 
Monterey  Institute  for  Research  in  Astronomy.  He  has  developed  a  stochastic 
Monte  Carlo  stellar  atmosphere  model  (SAM).  The  advantages  of  this  approach 
are  listed  below: 

1.  Interaction  of  photons  and  matter  can  be  modeled  explicitly  and  in  detail 
permitting  easy  inclusion  of  NLTE  effects; 

2.  Computations  are  done  in  spherical  geometry; 

3.  Influence  of  stellar  rotation  is  included; 
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4.  Influence  of  stellar  winds  is  included; 

5.  2  dimensional  geometiy  is  used  naturally  reproducing  limb  darkening;  and 

m 

6.  Output  is  in  2  dimensions;  X.  and  6  ;  where  6  is  the  an^e  from  the  surface 
normal. 

SAM  is  a  true  simulation,  with  the  fundamental  atomic  processes  being 
modeled.  The  advantage  of  this  approach  allows  the  modeler  to  investigate  funda¬ 
mental  physical  processes,  leaving  the  radiative  transfer  computations  to  the  com¬ 
puter. 

The  theoretical  concept  behind  SAM  is  simple.  At  large  optical  depth,  the 
radiation  field  is  acoirately  described  by  the  Planck  function.  Individual  photons 
selected  from  the  Planck  distribution  are  foUowed  as  they  migrate  through  the 
atmosphere.  The  properties  of  the  photosphere  are  then  accnirately  described  by 
the  LTE  assumptions  except,  as  stated  before,  for  the  radiation  field.  But  by  keep¬ 
ing  track  of  the  photons  firom  the  deeper  region  of  the  atmosphere,  where  the 
Planck  function  is  valid,  an  accurate  picture  of  the  radiation  field  can  be  produced. 
In  other  words,  SAM  is  able  to  reproduce  both  the  radiation  field  as  well  as  the 
state  of  the  atoms  in  the  important  LTE  region. 

If  one  then  continues  into  the  NLTE  region,  the  radiation  field  deter¬ 
mined  in  LTE  will  help  to  formulate  the  characteristics  of  the  NLTE  portion  of  the 
atmosphere,  thus  approximating  a  solution  to  the  NLTE  radiation  transfer  prob¬ 
lem. 

Stellar  wind  effects  can  also  be  simulated.  The  outer  atmosphere  of  SAM 
can  easily  be  set  in  motion.  This  allows  modelers  to  gain  insist  into  the  effects  of 
rapidly  moving  atmospheres  on  line  profiles. 
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2.  SAMI  13 


SAM1_13,  the  most  recent  version  of  SAM,  is  used  as  the  starting  point 
for  this  thesis.  The  method  employed  by  SAM1_13  to  simulate  the  radiative 

e 

transfer  problem  is  to  generate  a  ntunber  (usually  5(X),0(X)  photons/A)  of  photons 
and  distribute  these  photons  evenly  over  the  wavelengths  in  question.  The 
assumption  being  the  Planck  function  is  flat  over  a  few  angstroms. 

Photons  are  stalled  at  a  calculated  thermalization  depth,  usually  deep 
within  the  photosphere.  The  calculations  of  thermalization  depth  are  based  on  con* 
tinuiun  photons  with  optical  depth  of  10.  The  individual  photons  are  now  permitted 
to  migrate  through  the  atmosphere.  Random  absorption  probabilities  are  produced 
for  continuum  and  line  processes.  Since  continuous  variation  of  atmospheric  pro¬ 
perties  is  impractical,  atmospheric  parameters  are  held  constant  for  one  to  ten  per¬ 
cent  absorption  probability.  These  small  absorption  probabilities  are  added  until 
the  initial  randomly  chosen  absorption  probability  is  reached.  This  procedure  is 
followed  for  both  line  and  continuum  processes.  The  lengths  traveled  for  the  two 
processes  are  then  compared,  with  the  shortest  being  chosen  by  the  computer  as 
the  absorption  process  the  photon  imdergoes.  At  this  point  the  photon  is  con¬ 
sidered  absorbed  and  a  re-emitted  photon  is  produced.  For  simplicity,  coUisional 
processes  are  not  considered.  The  new  photon  is  redistributed  in  wavelength  and 
emission  angle  according  to  the  absorption  process  involved  (i.e.,  H~  absorption, 
line  absorption).  This  process  is  continued  for  the  photon  until  an  assumed  escape 
altitude  is  reached,  usually  the  top  of  the  photosphere,  or  it  is  destroyed  by  pro¬ 
pagating  back  into  the  thermalization  depth.  Radiated  photons  are  then  collected 
into  bins  chosen  by  the  user,  typically  of  width  .01  angstroms.  Angle  to  the  surface 
normal  at  the  time  of  departure  is  also  retained  and  radiated  photons  are  binned 
accordingly. 

SAM1_13  requires  two  programs  to  be  run  prior  to  actual  photon  genera¬ 
tion.  These  programs  read  in  atmospheric  parameters  and  generate  constants  that 
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are  needed  by  SAM.  Figure  9  ehows  the  order  in  which  the  different  programs  are 
run.  Figure  10  shows  the  conceptual  flow  diagram  of  SAM’s  logic.  Figure  11  is  a 
detailed  diagram  of  the  photon  absorption  loope  as  described  above. 


SAM  PROGRAM  SEQUENCE 


Figure  9:  Order  of  execution  of  programs. 
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SAM  CONCEPTUAL  FLOW  DIAGRAM 


Input  spectral  region  of  interest  and 
atomic  parameters 

Begin  calculation  cycle  (blackbody  pho- 

tons  +  collisional  photons) 

I 


Start  blackbody  photons  k 


a)  Calculate  propagation  distance  for 
line  and  continuum  processes _ 


b)  Does  photon  escape? 


yea 


no 


152  c)  Does  photon  re-thennaliae? 


ye» 


d)  Record  photon  X,9  K- 


Is  blackbody  cycle  complete? 


no 


yes 


Start  coUisional  photons  If 


Same  step  as  a-d  al»ve 


Is  collisional  cycle  complete? 


yes 


Start  new  calculation  cycle 


Figure  10:  The  logic  flow  of  SAM  during  radiative  transfer. 
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LINE  ABSORPTION  LOOP 


CONTINUUM  ABSORPTION  LOOP 


Figure  11:  The  logic  flow  of  the  two  absorption  loops  in  SAM. 


Spherical  geometry  is  the  natural  choice  for  photons  propagating  in  SAM. 
The  position  of  photon’s  can  be  described  by  the  three  spherical  cc  irdinates  p,  6,  , 
and  .  Stellar  atmospheres  have  spherical  symmetry  rendering  p  the  only  impor¬ 
tant  coordinate  in  describing  the  photon’s  position.  The  photon’s  direction  of  pro¬ 
pagation  can  be  described  by  the  two  coordinates  6  and  ^  .  Where  6  is  the  angle 
firom  the  surface  normal  (p),  and  p  is  the  remaining  angle  defining  the  orientation 
of  the  plane  formed  by  0  and  p.  The  geometry  is  simplified  by  assuming  that  the 
orientation  of  this  plane  ip  angle)  is  symmetric  for  all  photons.  The  P  dependence 


19 


can  be  condensed  into  the  random  6  calculations.  With  the  condensation  of  the 
coordinate  system  becomes  a  floating  reference  firame. 


Geometzy  of  SAM 


Figure  12:  Orientation  of  p.  6,  and 

In  order  to  minimize  computation  time,  all  photons  being  radiated  are 
assumed  to  propagate  in  the  direction  of  the  (4>8erver.  At  first,  this  may  seem  to 
be  an  inappropriate  assumption,  however,  once  the  photons  have  left  the  atmo¬ 
sphere  of  the  star,  the  6  dependence  is  removed  and  photons  may  be  redistributed 
to  propagate  in  the  direction  of  the  observer.  Another  way  to  look  at  thic 
simplification  is  to  think  of  the  photons  traveling  backwards  finom  the  observer. 
When  the  photons  are  first  absorbed  at  the  upper  surface  of  the  star  the  incoming 
0  value  is  completely  lost.  It  is  at  this  point  where  0  can  have  any  value,  therefore 
a  0  is  chosen  so  the  photon  propagates  in  the  direction  of  the  observer. 
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3.  The  shortcomings  of  SAM1_13 

Although  SAM1_13  generates  some  interesting  restilts,  it  has  several 
shortcomings.  1)  It  fails  to  accoiint  for  an  unequal  number  of  photons  being  pro¬ 
duced  and  destroyed  in  the  atmosphere  above  the  thermalization  depth  by  coUi- 
sional  processes.  2)  It  fails  to  account  for  energy  conservation  in  the  upper  photo¬ 
sphere.  This  failure  caused  the  upper  atmosphere  atoms  to  have  LTE  population 
levels  when  NLTE  population  levels  would  be  more  accurate.  3)  It  DeuIs  to  have  a 
mechanism  for  imposing  stellar  rotation  effects  onto  the  radiated  photons.  Rota¬ 
tion  effects  can  be  significant  in  some  stars  (i.e.,  line  broadening  of  up  to  several 
angstroms).  And  4)  it  contains  simplified  wavelength  redistribution  functions.  In 
addition,  it  is  felt  that  due  to  the  large  quantity  of  photons  being  generated  and 
the  many  interactions  each  photon  imdergoes,  a  more  rapid  method  of  determining 
current  atmospheric  parameters  is  needed. 


21 


n.  OBJECTIVES  OF  TfflS  THESIS 


The  main  objective  of  this  thesis  is  to  improve  SAM  by  including  the  following 
physical  processes  to  provide  a  more  complete  model,  including  some  NLTE  effects. 
The  physical  processes  being  added  are:  stellar  rotation,  collisional  excitation,  coUi- 
sional  de-exdtation,  optimized  photon  starting  depth,  and  an  improved  photon 
redistribution  fimction.  Additionally,  minimizing  the  time  required  to  calculate  stel¬ 
lar  atmosphere  parameters  will  also  be  added.  Once  the  improvements  are  made, 
SAM-generated  lines  are  tested  against  tiiree  observed  solar  lines.  These  lines  are 
Ti  X6258.1  A,  3d^4s  -  ZdHp ,  multiplet  #  104,  J  4-5;  K  X4044.1  A,  48-5p,  multiplet 

#  3,  J  -  1-^;  Sr  X4607.3  A,  5s*  -5s5p,  multiplet  #  2,  J  0-1.  In  addition,  a  com- 

paiison  of  the  SAM-generated  limb-darkening  with  the  observed  solar  limb- 
darkening  is  made.  Finally,  an  attempt  is  made  to  reproduce  the  characteristic 
structure  observed  in  the  spectra  of  stars  with  strong  non-thermal  winds. 
Although  NLTE  features  are  programmed  into  SAM  (i.e.,  two  temperatures  for  the 
outer  photosphere),  no  attempt  is  made  in  this  work  to  generate  NLTE  outer  atmo¬ 
sphere  atomic  population  levels. 
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in.  IMPROVEMENTS  TO  SIMULATION 


A.  STELLAR  ROTATION 

Doppler  shifting  due  to  stellar  rotation  can  significantly  broaden  observed  line 
profiles.  Since  this  is  due  to  the  integration  of  shifted,  unbroadened  lines  fttun  all 
portions  of  the  disk  visible  to  the  observer,  the  total  amount  of  lig^t  removed  finm 
the  continuum  is  unchanged.  At  high  rotational  velocities,  a  few  himdred  kilome¬ 
ters  per  second,  additional  line  broadening  can  be  as  much  as  several  angstroms. 
Inclusion  in  SAM  is  necessary  to  corrsctly  reproduce  line  profiles  fi;t>m  stars  rotat¬ 
ing  faster  then  1  km/sec.  The  star  is  assumed  to  rotate  as  a  solid  body.  This 
assiimption  is  justified  for  many  stars.  For  example  the  difierential  rotation  for 
the  sun  from  the  equator  to  the  75"  latitude  is  only  25%.  Because  SAM  follows 
individual  photons  as  they  travel  through  the  rotating  stellar  atmosphere,  Doppler 
shifts  within  the  atmosphere  (internal  Doppler  shifts)  must  be  examined  as  well  as 
those  affecting  the  observations  of  an  external  observer  (external  Doppler  shifts). 

There  are  no  internal  Doppler  shifts  because  of  the  solid  body  assxunption. 
Since  two  internal  points  in  the  atmosphere  have  no  relative  motion,  no  Doppler 
shifting  can  occur  to  a  photon  traveling  between  these  two  internal  points. 

The  only  broadening  due  to  rotation  comes  firom  external  Doppler  shifts.  Fig- 
xue  13  shows  photons  escaping  fi:x)m  the  limbs  of  the  star  which  will  appear  to  the 
observer  as  blue  and  red  shifted  with  respect  to  the  center  of  the  disk. 

Each  exiting  photon  in  SAM  must  be  Doppler  shifted  according  to  its  exit 
location.  The  output  of  SAM  is  divided  into  angular  increments.  A  rotating  body 
inclined  at  an  angle  i  to  the  line  of  sight  of  the  observer  can  have  its  surface 
divided  into  a  series  of  strips  of  imique  radial  velocities.  Superimposed  over  the 
angular  divisions  of  SAM  are  these  unique  radial  velocity  strips.  In  a  separate 
software  routine  called  VEL,  the  geometrical  regions  are  calculated  and  the  indivi¬ 
dual  photons  firom  each  region  are  then  Doppler  shifted  as  appropriate.  The 
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Photon  Doppler  shift  diagram 


Figure  13:  Photons  are  Doppler  shifted  as  seen  by  the  observer  depending  on 
the  portion  of  the  star  from  which  they  radiate. 


general  equation  for  this  procedure  is 


X 

l+v^(v)  ’ 


(4) 


where  NPCv^r)  =  number  of  photons  as  a  function  of  v  and  r, 

Ap  =  ratio  of  total  area  to  partial  area  of  rings  and  velocity  strips,  and 
Vf^  =  ratio  of  strip  velocity  to  speed  of  light 


The  input  rotational  speed  for  VEL  is  equal  to  v  sin(i)  .  Figure  14  shows  the 
geometric  regions  of  interest  and  the  complicated  areas  that  must  be  calciilated. 
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Velocity  and  angular  regioxis 


Figure  14:  Velocity  and  angular  regions  as  projected  onto  the  surface  of  the 
star.  The  regions  that  must  be  calculated  are  defined  by  the  intersection  of 
the  these  geometric  regions. 

O 

Figure  15  shows  the  resiilts  of  applying  VEL  to  Ti  X6258.1  A  for  three  different 
rotational  velocities.  These  resxUts  agree  with  theory  and  show  that  rotation 
effects  can  be  quite  significant.  Solar  rotation  is  2  km/sec,  and  has  little  effect  on 
line  profiles. 
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Stellar  rotational  broadenine 


Figure  15:  Different  rotational  velocities  as  applied  to  Ti  X6258.1A  The  top 
line  is  for  a  rotational  velocity  of  50  km/sec,  the  middle  line  is  for  rotational 
velocity  of  10  km/sec,  and  the  bottom  line  is  for  1.0  km/sec.  The  flux  ofi&et  is 
induced  for  viewing  clarity. 

B.  COLUSIONAL  EXCITATION  AND  DE-EXCITATION 

A  msgor  step  towards  making  SAM  a  tiilly  NLTE  simulation  is  the  introduc¬ 
tion  of  colUsional  processes.  Photons  need  to  be  added  as  a  result  of  coUisional 
excitation  of  an  atom  followed  by  photo-de-excitation.  In  addition,  photons  need  to 
be  removed  as  a  result  of  photo-exdted  atoms  which  iindergo  coUisional  de- 
exdtation. 

The  key  to  each  of  these  processes  is  the  value  of  the  coUisional  exdtation 
(C^)  and  de-exdtation  (Cui)  rates.  From  Jefferies  [Ref.  16],  a  semi-empirical 
equation  based  on  the  dipole  approximation  provides  a  means  to  calculate  the  exd¬ 
tation  rate  (Qu )  from  the  lower  to  the  upper  state. 
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Ciu  =2.16a-^“exp(-<x)T-iVne 


(5) 


where  a  =  -==  , 

KT 

T|o  =  energy  of  excitation, 
rig  =  electron  density, 
f  =  oscillator  strength,  and 
T  =  temperature. 

With  the  value  of  Ciu  known,  calculation  of  C^i  follows  from  the  principle  of 
detailed  balance.  This  principle  states  that,  in  TE,  all  processes  occur  exactly  as 
rapidly  as  their  inverse.  Although  TE  is  required  to  derive  this  equation  (eq  6),  the 
fact  that  it  only  contains  various  atomic  parameters  that  are  properties  only  of  the 
atoms  themselves  allows  the  equation  to  be  valid  under  NLTE  conditions  as  well 
The  equation 

(6) 

8u  ^ 

where  gi  «  statistical  weight  lower  level,  and 
gu  -  statistical  weight  upper  level 

can  be  used  to  relate  these  two  coUisional  rates. 

To  compute  the  niimber  and  starting  depth  of  the  photons  produced  by  the 
coUisional  excitation  process,  the  sum  of  all  the  photons  produced  by  this  coUisional 
process  is  first  calculated  firom  the  starting  depth  of  the  blackbody  photons  (xq)  to 
the  exit  altitude  of  the  atmosphere  (Zf ).  The  equation  for  the  total  number  of  coUi¬ 
sional  photons  from  depth  Xq  to  x/-  is 

(7) 

*0 
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and  for  the  blackbody  photons  is 


2ffc 


e^-1 


(8) 


where:  Nc^  =  Total  niunber  of  collisional  photon  /cm^  firom  ZotoZf  ,  and 
Nfcfc  =  Number  of  blackbody  photon  !  cm}  . 

In  this  way  a  ratio  of  the  number  of  blackbody  photons  to  coUisionally- 
produced  photons  can  be  found.  For  example,  if  the  ratio  of  blackbody  photons  to 
ooUisionally-produced  photons  is  1  to  3,  for  every  blackbody  photon  SAM  generates, 
it  must  also  produce  three  collisionally-produced  photons.  A  collisional  photon  vs. 
depth  distribution  is  also  calculated  during  the  siunming  of  the  collisional  photons. 
Collisionally-produced  photons  are  assigned  a  starting  depth  randomly  based  on 
this  probability  distribution. 

Photon  removal  caused  by  collisional  de-exdtation  is  calculated  in  the  line 
interaction  subroutine.  If  a  photon  is  absorbed  by  an  atom,  then  a  collisional  de- 
exdtation  probability  is  checked  against  a  random  number  and  those  photons  that 
fall  within  the  de-exdtation  probability  are  removed. 

In  this  way,  the  collisional  source  function,  dependent  on  depth,  and  its 
accompanying  extinction  process  are  simxilated.  Unlike  typical  treatments  by  radi¬ 
ative  transfer  theory,  SAM’s  source  function  does  not  need  to  include  photo- 
exdtation  and  de-exdtation.  These  processes  are  inherent  by  the  way  SAM  follows 
individual  photons. 

Figures  16  and  17  show  the  SAM-generated  Ti  X6258.1  A  line  with  and 
without  the  collisional  source  function  induded.  Note  the  line  without  the  colli¬ 
sional  source  function  (Figure  16)  is  deeper  than  the  line  with  the  source  function 
induded.  This  effect  is  caused  by  an  imbalance  between  photons  being  produced 
by  collisional  exdtation  and  photons  being  destroyed  by  collisional  de-exdtation.  In 
the  outer  photosphere  there  are  more  photons  added  by  collisional  exdtations  than 
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removed  by  ooUisional  de-exdtation. 

Ti  line  with  no  collisionally  excited  photon 


Figure  16:  SAM  generated  line  with  no  oollisionally-produced  photons  added. 
Ti  line  with  collisionally  excited  photons 


Figure  17:  SAM  generated  line  with  coUisionally-produced  photons  added. 
Note  that  the  line  depth  has  decreased  due  to  the  increased  number  of  photon 
finom  coUisional  excitation. 
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C.  INTERPOLATION  TABLES 


In  an  effort  to  reduce  the  program’s  computation  time,  an  analysis  of  SAM’s 
computation  properties  was  xindertaken.  The  results  indicated  the  computer,  a 
S\m  SPRACl  workstation,  was  spending  a  significant  amount  of  CPU  time  calcu¬ 
lating  and  re-calculating  atmospheric  parameters.  Specifically,  the  subroutine 
STAR,  where  spline  fits  to  atmospheric  parameters  were  being  calculated,  was  the 
largest  user  of  CPU  time. 

A  simple  re-design  of  the  subroutine,  replacing  the  spline  fits  with  interpola¬ 
tion  tables,  was  completed.  The  interpolation  tables  are  calculated  in  NSUM  using 
the  same  spline  fit  routine.  The  interpolation  tables  are  incremented  in  1  km  dis¬ 
tances.  This  distance  may  be  varied  when  stars  other  then  the  sun  are  examined. 
To  ensure  no  detail  in  the  structure  of  the  atmosphere  is  lost,  graphs  of  all  atmos¬ 
pheric  parameter  were  compared  to  the  interpolation  tables  (see  Appendix  A).  The 
net  effect  of  this  re-design  is  a  300%  increased  in  algorithm  efficiency  with  negligi¬ 
ble  decrease  in  acciuracy. 

D.  OPTIMIZATION  OF  PHOTON  STARTING  DEPTH 

Another  significant  degradation  to  computation  speed  is  the  chosen  photon 
starting  depth.  It  was  found  that  photons  which  started  too  deep  within  the  star 
would  xmdergoes  needless  absorption  and  re-emission.  A  starting  depth  which 
minimized  computer  time  is  in  order.  The  depth  chosen  is  where  the  line  absorp¬ 
tion  and  continuum  absorption  processes  had  identical  optical  path  lengths.  Above 
this  point  in  the  atmosphere,  the  difference  between  these  two  optical  path  lengths 
shaped  the  outgoing  photon  flux.  Below  this  point,  the  differences  in  the  path 
length  may  shape  the  outgoing  flux,  but  the  flux  is  completely  re-ordered  as  it 

O 

migrates  through  the  equal  path  length  region.  Figure  18  shows  Ti  (X6258.1  A  ) 
line  generated  at  the  optimal  depth  (484km)  and  one  generated  with  z  equal  to 
490km. 
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There  is  little  different  between  the  two  line  profiles,  but  the  computer  run 
for  the  line  profile  with  z  equal  to  484km  is  10%  as  fast 


Ti  line  generated  at  two  different  starting  depths 


Figure  18:  The  same  Ti  line  started  at  two  different  starting  depth.  Note  that 
the  line  shape  is  almost  identical.  Difference  in  fiux  of  the  two  lines  is  in¬ 
duced  for  viewing  clearity. 

E.  REDISTRIBUTION 

Photons  undergo  many  different  absorption  and  scattering  processes  while 
propagating  through  the  stellar  atmosphere.  How  these  photons  are  redistributed 
in  fiiuquency  is  of  great  importance  in  determining  the  final  line  profile.  SAM1_13 
had  two  simplified  redistribution  functions;  complete  Compton-scattering  redistri¬ 
bution  for  the  continuum  absorption  process,  and  modified  complete  Doppler  redis¬ 
tribution  for  the  line  absorption  process.  These  two  simplification  are  unsuitable 
for  accurate  line  profile  simulation.  Improvement  to  both  finequency  redistribution 
functions  is  necessary. 
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The  migor  source  of  continuum  absorption  in  the  sim  is  firom  H~  .  H~  has  a 
low  ionization  energy  at  .754  ev,  therefore  photons  with  wavelengths  less  then 

O 

16,450  A  have  enough  energy  to  ionize  H~  .  Compared  to  other  sources  in  the  visi¬ 
ble  region  of  the  continuum,  H~  is  stronger  by  a  factor  of  4  than  that  due  to  neu¬ 
tral  H  and  much  larger  than  that  for  other  processes.  Therefore,  for  the  sun,  all 
continuum  absorption  is  assumed  to  be  caused  by  H~  .  This  assumption  is  valid 
for  many  types  of  stars.  H~  low  ionization  energy  results  in  a  simple  redistribution 
function.  All  photons  absorbed  by  the  continuum  process  are  completely  redistri¬ 
buted  in  firequency.  In  consideration  of  minitniring  computation  efforts,  photons 
redistributed  out  of  the  wavelength  region  of  interest  are  replaced  by  the  introduc¬ 
tion  of  a  random  wavelength  photon  firom  these  other  regions.  In  thia  way,  photon 
flux  is  conserved. 

The  line-absorption  frequency-redistribution  fimetion  is  a  combination  of 
Doppler  redistribution  (Gatissian  in  shape)  and  coUisional  redistribution 
(Lorentzian  in  shape).  It  is  the  combination  of  these  two  distributions  that  result 
in  the  familiar  Voigt  profile  for  line  emission  and  absorption.  The  rate  of  collisions, 
the  Voigt  a  parameter,  the  kinetic  temperatiire,  and  the  spontaneous  emission 
coefficient  (Einstein  A),  are  all  important  factors  in  determining  line-absorption 
firequency-redistribution.  Doppler  redistribution  is  based  on  atoms  having  Maxwel¬ 
lian  velocity  distribution  [Ref  17].  The  standard  Gaussian  profile  is  given  by 

p-**/ao* 

=  (9) 

o 

where:  o  =  half  width  at  half  maximum. 


The  Doppler  redistribution  function  from  Hummer  [Ref  18]  is  given  by 


4(v)  = 


If  we  let  Af  =  then 


(10) 
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(11) 


where  »  Aj  *  ’ 

Xo  =  initial  wavelength, 
c  =  speed  of  U^t, 

K  s  Boltzman  constant, 

T  s  kinetic  temperature, 

M  s  mass  of  the  atom, 
s  average  Doppler  shift. 

For  the  collisional  redistribution,  the  standard  Lorentzian  profile  is  given  by 

P(x)= - (12) 


The  collisional  redistribution  function  from  Hummer  [Ref.  18]  is  given  by 

^v)  = _ - 

x((Av)*4(A’‘+A*)*) 

where  A"  =  average  natural  broadening, 

A'  =  average  collisional  broadening. 

Let  o  =  A^+A®  and  rearranging 

1 


(13) 


<t(v)  = 


JCO(l+^) 


=  P(x). 


(14) 


The  Voigt  a  parameter  is 


a  - 


therefore  a^A^a  . 
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Normal  random  niunber  distribution  was  modified  by  the  Boz-Muller  method 
[Ref-  19]  to  produce  a  temperature-dependent  Gaxissian  (eq  15).  Random 
Lorentzian  distribution  dependent  on  the  Voigt  a  parameter  and  temperatm*e  was 
produced  by  inverse  transform  method  (eq  16)[Ref.  19]. 

X  shift  =  V-2/n  (r  )cos  (2nr  yJ  —  (15) 

^  M  c 

X  shift  -  tan  (nr  ^  a  (16) 

where  r  =  random  number  between  zero  and  one. 

These  redistribution  functions  were  tested  by  generating  a  million  random 
n\imbers  and  plotting  the  resulting  X  shift  (Figure  19  and  20). 

Gaussian  redistribution  function 


1100 .00 


Figure  19:  Gaussian  curve  produced  by  redistribution  function. 
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Lorentzian  redistribution  function 


4100.00 


Figure  20:  Lorentzian  ciuve  produced  by  redistribution  ftmction. 

Since  collisional  redistribution  is  dependent  on  whether  an  atom  undergoes  a 
collisional  interaction,  not  all  line-absorbed  photons  are  collisionally  redistributed. 
As  a  first  order  approximation,  it  can  be  assumed  that  the  ratio  of  collisional  de¬ 
excitation  compared  to  total  de-excitation  is  the  same  as  the  ratio  of  photons  which 
undergoing  collisional  redistribution  compared  to  just  Doppler  redistribution. 

It  is  with  these  two  redistribution  fimctions  (eq  15  and  16)  that  re-emitted 
photons  are  redistributed  in  wavelength  in  SAM.  The  photons,  once  absorbed,  are 
re-emitted  with  a  new  wavelength  assigned  randomly  by  a  combination  of  these 
two  functions.  These  functions  are  a  mqjor  improvement  over  the  simplified  func¬ 
tions  being  used  in  SAM1_13. 
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IV.  TESTING  AND  RESULTS 


Three  methods  of  testing  the  validity  of  SAM  were  undertaken.  The  first 
method  involved  comparisons  of  theoretical  mean  firee  path  lengths  to  computed 
mean  firee  path  lengths  in  a  homogeneoxxs  atmosphere.  The  second  method  is  com¬ 
parison  of  observed  absorption  line  profiles  to  SAM-generated  profiles.  Finally, 
observational  limb  darkening  is  compared  to  SAM-generated  limb  darkening.  To 
demonstrate  the  fiexibility  of  SAM,  stellar  turbulence  and  stellar  wind  effects  are 
applied  to  the  solar  Ti  X6258.1  A  line. 

A.  MEAN  FREE  PATH 

In  this  test  both  the  line  absorption  subroutine  and  the  continuum  absorption 
subroutines  are  tested.  Although  stellar  atmospheres  are  not  homogeneous,  the 
homogeneous  assumption  allows  mean  firee  path  calculations  to  be  made.  There¬ 
fore  this  is  an  important  test  for  SAM.  Mean  free  path  calculations  arc  shown 
below. 

Treating  the  atmosphere  as  homogeneous  and  including  only  absorption 
processes  of  the  continutun,  the  expected  normal  exponential  decay  equation  is 
obtained. 

^  =  exp  (-Kvplr )  (17) 

*0 

where:  I  =  intensity  at  distance  L  through  the  medium 
7o  =  initial  intensity, 

Ky  =  continuum  opacity  per  gram,  and 
p  =  density. 

The  average  length  U>aoe )  is; 


where  and  p  are  given  by  the  photospheric  model  being  used. 

For  the  mean  free  path  of  the  line  absorption  process,  the  equation  is; 


Therefore, 


where 


^=e*p(-iCvL). 


(19) 


(20) 


gufNjicr^ 

me 


HiuA) 


f  =  oscillator  strength, 

gu  =  statistical  weight  upper  level, 

U  -  partition  function  of  the  atom, 

Ni  =  population  level  of  i  state, 

0  =  5039.99/kinetic  temperature, 

X  =  excitation  energy  of  lower  level  (eV), 

Avx)  =  Doppler  shift,  and 
H{u  /I )  =  Hjerting  function. 

SAM  generated  values  were  produced  by  averaging  3000  path  lengths 
through  the  homogeneous  medium.  Table  1  and  2  contains  various  input  parame¬ 
ters  as  well  as  both  calculated  and  SAM-generated  .  Mean  free  path 
comparisons  for  the  continuum  subroutine  are  done  for  four  different  continuum 
optical  depths,  while  the  line  absorption  subroutine  are  done  for  four  different  Une 
optical  depths. 
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TABLE  1  rCONTINUUM 


X 

depth  (km) 

' 

ix-S-) 

cm® 

LoMikm  Xtheory] 

Loue(^XSAM] 

.001 

83 

1.28x10-2 

1.37x10-® 

56800 

56300 

.01 

213 

3.61jc10-2 

4.88x10-® 

5670 

5630 

.1 

344 

l.lOtrlO-^ 

1.59X10-2 

569 

565 

1.0 

455 

.945 

3.31x10-2 

31.9 

33.9 

m 

500 

17.5 

3.15x10-2 

1.81 

1.90 

TABLE  2  :LINE 

B 

depth  (km) 

H(u.a) 

Kv/cm-^) 

L««(^  Xtheory] 

L«^(ikmXSAM] 

.05 

55 

.9628 

9.09x10-® 

1100 

1109 

.18 

83 

.9497 

2.14x  10-® 

467 

463 

3.17 

213 

.8569 

1.49x10-2 

67.2 

66.0 

15.0 

344 

.6079 

4.35x10-2 

22.9 

22.9 

A  comparison  of  the  results  reveals  the  good  agreement  between  the  typical 
mean  free  paths  calculated  by  SAM  and  those  predicted  by  theory.  The  minor 
difference  between  the  two  values  may  by  accounted  for  by  the  averaging  of  the 
random  processes. 
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B.  COMPARISON  TO  OBSERVATION 


The  second  method  of  validation  is  the  comparison  of  the  computer  generated 
lines  to  observed  lines.  Line  observations  are  taken  from  the  National  Solar  Obser¬ 
vatory  Atlas  [Ref.  20]. 

1.  Convergence 

Figures  21,  22  and  23  show  a  time  series  firom  SAM.  Since  SAM  is  a 
Monte  Carlo  simulation,  the  more  photons  SAM  generates,  the  higher  the  resiilting 
signal  to  noise  ratio  in  the  SAM  generated  curve  becomes.  The  trade  off  can  be 
seen  in  these  three  figures.  All  SAM  line  profiles  were  generated  on  a  Sun 
SPRACl  workstation  and  CPU  times  stated  are  from  this  computer.  Figure  21 
required  only  2  minutes  of  CPU  time  and  produced  10,000  photons.  Figure  22 
required  20  minutes  of  CPU  time  and  produced  100,000  photons.  Figure  23 
required  about  200  minutes  of  CPU  time  producing  1,000,000  photons. 


10,000  photons 


wavelength 

Figure  21: 10,000  photons  produced  by  SAM  in  2  xninutes 
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100,000  photons 


Figiire  22: 100,000  photons  produced  by  SAM  in  20  minutes. 

1,000,000  photons 


Figure  23:  1,000,000  photons  produced  by  SAM  in  approximately  200 


minutes. 

Standard  fluctuation  calculation  reveals  that  the  fluctuation  per  bin  is  in 
agreement  with  random  processes. 


F  =^INavt 


(20) 
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where:  F  =  Fluctuation  per  bin,  and 

^ave  -  Average  niunber  of  photons  per  bin. 

2.  Lone  ssmthesis 

One  of  the  m^or  applications  of  SAM  is  the  fitting  of  stellar  and  atomic 
parameters.  For  this  thesis,  three  solar  line  are  synthesized.  The  portions  of  the 
solar  spectrum  being  synthesized  are  shown  in  Figiures  24,  25  and  26.  The  figures 
show  two  sets  of  lines.  The  dark  line  is  the  intensity  of  the  solar  flux  on  a  scale 
between  0  and  100%.  The  lighter  line  is  the  more  detailed  scale  running  firom  90 
to  100%.  It  is  the  darker  line  that  SAM  is  modeling.  The  wavelengths  given  in 
these  figures  are  in  nanometers.  These  lines  were  chosen  principally  because  of 
the  availability  of  the  two  atomic  constants;  spontaneous  emission  coefficient  (Ein¬ 
stein  A),  and  quadratic  stark  interaction  constant  (C4).  For  comparison  to  SAM 
generated  lines,  the  absorption  lines  near  the  absorption  line  being  modeled  have 
been  replaced  by  a  smoothed  continuum. 

Observed  portions  of  solar  spectrum 


Figure  24;  Observed  solar  spectrum  near  Sr  X4607.3  absorption  line. 
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404.3  404.4 

Figure  25:  Observed  solar  spectrum  near  K  X4044.1  absorption  line. 


Figure  26:  Observed  solar  spectrum  near  Ti  X6258.1  absorption  line. 

Initial  comparisons  to  estimated  abundance,  and  C4  values  are  shown  in 
Figures  27,  28,  and  29.  Notice  that  most  of  the  fits  lack  the  required  line  strength 
and  line  broadening.  Abundance  values  are  given  as  the  log  of  the  ratio  of  atom 
abundance  to  hydrogen  abundance,  with  log  H  abundance  =  12.0. 
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flux 


Comparison  of  observed  lines  to  SAM-generated  lines 
using  initial  C4  and  abundance  values 


Figure  27.  Observed  K  X4044.1  line  to  SAM-generated  line  with  abun¬ 
dance  *  5.6  and  log  C4  =  -13.59. 


Figure  28:  Observed  Sr  X  4607.3  line  to  SAM-generated  line  with  abun¬ 
dance  =  2.7  and  log  C4  =  -15.81. 
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Figure  29:  Observed  Ti  X  6258.1  line  to  SAM-generated  line  with  abun¬ 
dance  =  4.78  and  log  C4  =  -10.41. 

The  initial  values  of  abundance  and  C4  are  adjusted  until  a  good  fit  is 
obtained.  Adjustment  is  appropriate  as  these  values  typically  have  a  large  uncer¬ 
tainties.  Figures  30,  31,  and  32  show  SAM-generated  fits  with  adjusted  atomic 
values.  Quality  of  the  fits  can  be  compared  to  that  in  Figure  3. 
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flux  flux 


Comparison  of  observed  lines  to  SAM-generated  lines 
using  a4justed  C4  and  abundance  values 


Figure  30:  Observed  K  X4044.1  line  to  SAM-generated  line  with  abiin- 
dance  =  5.81  and  log  C4  =  -11.19. 


Figure  31:  Observed  Sr  X4607.3  line  to  SAM-generated  line  with  abun¬ 
dance  =  3.00  and  log  C4  =  -9.81. 
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Figure  32:  Observed  Ti  X6258.1  line  to  SAM-generated  line  with  abun¬ 
dance  =  4.80  and  log  C4  =  -10.41. 

C.  LIMB  DARKENING 

Solar  radiation  intensity  ifirgslcm^)  decreases  firom  the  center  of  the  solar 
disk  to  the  limb.  This  limb  darkening  effect  is  caused  by  successive  viewing  of 
higher  and  higher  photospheric  layers,  where  temperatinre  is  lower  and  the  soxirce 
function  is  less  intense.  It  is  the  limb  darkening  effect  that  allows  temperature 
distributions  and  approximate  source  functions  to  be  determined  as  a  function  of 
height  in  the  stellar  atmosphere.  These  are  important  parameters  for  stellar 
modeling.  The  imderstanding  and  successful  simulation  of  this  effect  is  therefore 
apphcable  and  relevant. 

The  SAM  simulation  retains  two  pieces  of  information  on  exiting  photons: 
photon  wavelength  and  exit  angle  from  the  surface  normal.  If  SAM’s  output  is  con¬ 
sidered  spread  uniformly  over  the  solar  disk  then  the  angular  binning  of  the  exit¬ 
ing  photons  gives  the  exit  location  of  the  photons.  But  limb  darkening  measure¬ 
ments  are  taken  from  a  thin  spatial  region  from  the  center  of  the  solar  disk  to  its 
edge.  These  measurements  in  effect  reduce  the  sim’s  intensity  into  a  2 
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dimensional  cross  section,  where  the  intensity  is  a  function  of  angular  position. 
But  observations  of  limb  darkening  are  made  in  terms  of  intensity  (photons lcm\ 
therefore  each  angular  region’s  exiting  photon  flux  must  be  divided  by  the  projected 
siuface  area  normal  to  the  observer  to  obtain  the  necessary  intensity.  Figure  33 

O 

shows  a  comparison  of  SAM-generated  limb  darkening  for  5000A  to  observed 
curve.  Agreement  is  good  until  the  SAM-generated  limb  intensity  becomes  greater 
than  observations  aroxind  70  degrees. 


SAM  limb  darkening  compared  to  observations 


sam  generated  limb  dorkening  ol  5000A 
Observed  limb  dorkening  ot  5000A 


ongle 

o  o 

Figure  33:  SAM  generated  limb  darkening  at  5000A  compared  to  5000A  ob¬ 
served  curve.  Deviation  is  still  under  investigation. 


Investigations  into  this  effect  seem  to  indicate  it  is  related  to  the  continuum 
opacity  iCy  but  earlier  mean  free  path  calculations  indicate  Xy  is  calculated  correctly. 


D.  STELLAR  TURBULENCE 


Turbulence  is  the  motion  of  the  atmospheric  gases  on  a  scale  that  is  large 
compared  to  atomic  dimensions  but  small  compared  to  the  size  of  the  star.  Eddies, 
convection  cells,  and  prominence-like  structures  are  all  forms  of  turbulence.  These 
motions  cause  additional  Doppler  shifting  to  the  photons  and  can  greatly  affect  the 
line  shape. 

The  velocity  Doppler  shifting  due  to  stellar  tiirhulence  occurs  in  SAM  during 
the  calciilation  of  the  line-absorption  coefficient.  SAM  does  not  currently  Doppler 
shift  the  photons  due  to  stellar  turbulence  during  re-emission.  This  would  be  fairly 
easy  to  correct  if  the  natiire  of  the  turbulence  is  known.  For  the  piirpose  of  this 
thesis  the  re-emission  correction  is  ignored.  The  broadening  caused  by  increasing 

O 

the  stellar  turbulence  velocity  to  the  Sr  X4607.3  A  is  shown  in  Figure  35. 

SAM  generated  stellar  turbulence 


Figure  35;  SAM  generated  stellar  turbulence  for  three  different  velocities. 
The  top  line  is  for  a  turbulence  velocity  of  6.0  km/sec,  the  middle  line  is  for 
4.0  km/sec,  and  the  bottom  line  is  for  2.0  km/sec.  Offsets  are  induced  for 
viewing  clarity. 
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E.  STELLAR  WINDS 


Stellar  winds  differ  from  turbulence  in  that  the  gases  driven  by  convection  or 
radiation  obtain  the  necessary  velocity  to  escape  from  the  star.  The  amount  of 
escaping  matter  can  be  quite  large,  typically  1  earth  mass  per  year.  The  effects  of 
this  escaping  matter  on  stellar  lines  is  not  well  known.  Although  some  general 
aspects  of  stellar  wind  effects  on  line  profiles  can  be  anticipated.  The  rapidly  mov¬ 
ing  matter  will  most  certainly  Doppler  shift  emitted  photons.  The  Doppler  shifting 
will  only  be  in  blue  direction  as  the  material  which  is  moving  away  finom  the 
observer  will  be  blocked  by  the  star.  Consequently,  line  profiles  shoxild  contain 
increased  flux  in  the  shorter  wavelengths.  Many  additional  effects  and  the  driving 
forces  of  these  winds  are  not  well  known.  This  is  one  of  the  motivations  behind  the 
creation  of  SAM. 

In  an  effort  to  test  the  stellar  wind  routine  in  SAM,  a  10  km/sec  maximum 
stellar  wind  was  imposed  on  the  sun.  Although  this  is  clearly  unrealistic,  the  goal 
was  to  reproduce  some  of  the  effects  of  stellar  winds  on  a  known  spectrum.  The 
wind  speed  varied  linearly  with  altitude,  reaching  a  maximum  of  10  km/sec  at  the 
top  of  the  photosphere,  the  equation  for  this  variation  is 

wind  velocity  =  — - —  maximum  wind  velocity  ,  (22) 

zmax 

where  z  =  distance  out  of  the  photosphere  frnm  reference  depth  (500km), 
zmax  =  distance  to  top  of  photosphere  from  reference  depth. 

Figure  36  shows  the  effect  of  this  wind  on  the  earlier  synthesized  Ti  absorption 
line.  The  large  emission  spike  at  the  shorter  wavelengths  is  the  same  effect 
observed  in  stars  with  large  stellar  winds  [Ref.  21].  This  is  an  exciting  initial 
resiilt.  To  emphasize  the  ease  in  which  SAM  was  modified,  only  four  additional 
lines  of  code  were  required  to  simulate  movement  of  the  outer  photosphere.  Two 
additional  wind  velocities  were  also  run  as  shown  in  Figures  37  and  38. 
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Stellar  wind  effects  on  absorption  lines 


Figure  36:  SAM  generated  absorption  line  with  10  km/sec  maximum  stellar 
wind. 


wovelength 


Figure  37:  SAM  generated  absorption  line  with  3  km/sec  maximum  stellar 
wind. 
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flux 


Figure  38:  SAM  generated  absorption  line  with  20  km/sec  maximum  stellar 
wind. 
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V.  CONCLUSION 


The  primary  goal  of  this  thesis  has  been  the  modification  and  testing  of  a  new 
photon-following  approach  to  stellar  radiative  transfer.  The  approadi  has  several 
msgor  advantages  over  standard  stellar  radiative  transfer  methods  and  has  only 
recmitly  been  possible  with  the  aid  of  fast  computers.  It  was  demonstrated  that 
this  technique  accurately  and  efficiently  reproduced  observed  solar  absorption  lines, 
including  the  effects  of  stellar  rotation,  stellar  turbulence,  coUisional  excitation,  and 
coUisional  de-excitation.  Stellar  winds  and  limb  darkening  were  reproduced 
although  problems  still  remain  in  these  areas.  Nevertheless,  the  successfiil  simvda- 
tion  of  these  effects  is  an  improvement  over  previous  methods.  One  of  the  major 
advantages  of  this  approach  to  stellar  radiative  transfer  is  many  vastly-difierent 
stellar  phenomena  can  be  modeled  and  incorporated  into  a  single  simulation. 

To  demonstrate  the  tisefulness  of  this  technique,  new  abundance  and  C4 
values  for  three  solar  lines  were  determined.  This  technique  can  also  be  used  to 
investigate  almost  all  phenomena  related  to  stellar  atmospheres  and  radiative 
transfer  (i.e.,  testing  various  model  atmospheres,  stellar  wind  theories,  stellar  tur¬ 
bulence,  etc.) 

The  results  described  in  this  thesis  used  fairly  well-behaved  stellar  absorption 
lines.  An  extension  of  this  thesis  would  be  to  generate  true  NLTE  atomic  popula¬ 
tion  levels  in  the  outer  photosphere.  If  this  were  accomplished,  the  complex  NLTE 
lines  could  be  simulated.  Some  additional  extension  area  are:  accurate  represen¬ 
tation  of  the  Hjerting  function  when  the  Voigt  a  parameter  exceeds  unity,  new 
redistribution  functions  that  separate  natural  and  coUisional  broadening,  and 
improved  velocity  profiles  for  both  steUar  turbulence  and  steUar  winds.  Although 
these  problems  are  non-trivial,  there  are  no  foreseeable  road  blocks  for  this  new 
technique  to  simulate  a  large  variety  of  complex  stellar  lines. 
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APPENDIX  A.  SOLAR  ATMOSPHERIC  PARAMETERS 


In  charpter  III  section  C  spline  fits  to  atmospheric  parameters  [Ref.  22]  were 
replaced  with  interpolation  tables.  Figures  39  through  42  show  a  comparison 
between  spline  fit  curves  and  interpolation  values  for  temperature,  density,  opacity 
and  gas  pressure  as  each  vary  with  depth.  Incremental  distance  of  5km  was 
chosen  to  allow  the  difference  between  the  two  curves  to  be  seen.  Better  agree¬ 
ment  is  actually  found  in  SAM  because  of  the  irse  of  1km  increments.  Note  an 
exception  to  the  5km  increments  is  Figure  41  where  an  incremental  distance  of 
2km  is  chosen  because  of  the  reduced  distances  being  plotted. 


Figure  39;  Temperature  curve  generated  by  spline  fit  compare  to  interpolation 
tables 
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Figure  41:  Opacity  ( 


g  ^ 


curve  generated  by  spline  fit  compare  to  interpola¬ 


tion  tables 
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